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ABSTRACT 

The Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS) has 
created a large homogeneous database of intermediate redshift (0.2 < z < 1.0) 
type la supernovae (SNela). The SNLS team has shown that correlations exist 
between SN la rates, properties, and host galaxy star formation rates. The 
SNLS SNIa database has now been combined with a photometric redshift galaxy 
catalog and an optical galaxy cluster catalog to investigate the possible influence 
of galaxy clustering on the SN la rate, over and above the expected effect due to 
the dependence of SFR on clustering through the morphology-density relation. 
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We identify three cluster SNe la, plus three additional possible cluster SNe la, and 
find the SN la rate per unit mass in clusters at intermediate redshifts is consistent 
with the rate per unit mass in field early-type galaxies and the SN la cluster rate 
from low redshift cluster targeted surveys. We also find the number of SNe la in 
cluster environments to be within a factor of two of expectations from the two 
component SNIa rate model. 

Subject headings: supernovae:general — galaxies:clusters:general 



INTRODUCTION 



The importance of understanding Type la supernovae (SNe la) has risen greatly in recent 
years because of their pivotal role in demonstrating the accelerated expansion of the un iverse 
flRiess et al. 19981 : IPerlmutter et al. 19991 : kstier et al. 20061 : IWood- Vasev et al. 2007^ SNe la 
are generally agreed to be carbon-oxygen white dwarfs which have accreted sufficient mass 
from their companion star to initiate a thermonuclear explosion. Despite this consensus, sev- 
eral models exist for the companion type, the accretion proce ss, the delay time between star 



formation and SN l a explosion, and the explosion mechanism (jHillebrandt fc Niemeyer 2000 



Hoflich et al. 20031 ). Photometric and spectroscopic observations, and explosion simul ations, 
are making inroads towards a coherent picture of SNIa events (jMazzali et al. 20071 ). An- 
other (albeit indirect) path of investigation is studying the SNIa rate in different stellar 
populations. 

Compilation of the Cappellaro et al. (1999) supernova catalog with infrared data led 
Mannucci et al. (2005) to conclude that SNII, SNIb/c, and SNIa rates per unit stellar mass 
are directly correlated with host morphology and (B-K) colour, and to infer a correlation with 
host star formation activity. Based on this, Scannapieco & Bildsten (2005) expressed the 
SN la rate per unit stellar mass as the sum of a "delayed" component from old and a "prompt" 
component from young stellar populations. They parametrized the two components as A and 
B, proportional to the mass and star formation rate (SFR) of a galaxy respectively, and also 
found the "prompt" component to account for the the discrepancy between low observed rates 
of cluster SNe la and the high cluster iron abundance (IMaoz fc Gal- Yam 20041 ). This "A+B" 
two componen t model has since be en confirmed at intermediate redshifts with the large 
SNLS catalog (jSullivan et al. 2006al ); the observations can also be matched with progenitor 
populations that have distributions of delay times as described by Mannucci et al. (2006) 
and Pritchet et al. (2007). 



From the two component model, the SN la rate in galaxy clusters is expected to be lower 
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than in the field due to the morphology-density relation (jPostman fc Geller 1984 ): cluster 
galaxies are predominantly of early-type with little or no star formation. However, since 
clusters contain only a small fraction of the stellar mass of the Universe, it is conceivable 
that some hitherto undetected influence in such exotic environments could affect the SNIa 
rate. For example, the fraction of binary stars, or the rate of mass accretion onto the white 
dwarf, could be enhanced. The recent disco yery of an enhanced nova rate in the core of 
elliptical galaxy M87 is evidence for the latter ([Madrid et al. 20071 ). A second example is th e 
detected SNIa rate enhancement in radio- loud early-type galaxies (IDella Valle et al. 20051 ); 
such galaxies tend to be very luminous ellipticals in the centers of large clusters. 

The Wise Observatory Optical Transient Survey (WOOTS) targeting 140 Abell clus- 
ters confirmed the SNe la rate per unit m ass in low redshift galaxy clusters to be consis- 
tent with the rate in early- type galaxies (ISharon et al. 20071 ). Recently, Mannucci et al. 
(2007) analyzed the Cappellaro et al. (1999) sample of 136 low redshift SNe (z < 0.04) 
and found the SNIa rate in cluster early-type galaxies is enhanced by a factor of > 3 
over field early-type galaxies. The large database of intermediate to high redshift SNe la 
compiled by the Ca nada-France Hawaii Telescope Supernova Legacy Survey (CFHT SNLS) 
( lAstier et al. 20061 ). an d the publicly ava il able photometric redshift galaxy and cluster cat- 
alogs for SNLS fields (lllbert et al. 20061 ; lOlsen et al. 20071 ) are ideal for extending these 
investigations to higher redshifts. As SNLS is not a cluster-targeted survey, further advan- 
tages over Sharon et al. (2007) include using a flexible parametrization of galaxy clustering 
in the environments of SNe la, and determining simultaneous field SN la rates for comparison 
to cluster rates. 

§ [2] describes the SNIa, galaxy, and cluster catalogs used in this experiment. §[3] and §|4] 
present two independent approaches to statistically evaluate the effects of clustering on the 
SNIa rate per unit mass. §|3]uses a continuous clustering strength estimator to compare the 
SNIa rate per unit mass inside and outside clustered environments. § H] identifies six SNIa 
in Olsen et al. (2007) galaxy clusters, considers the probability of these observations based 
on expectations from the two component rate model, and calculates the SN la rate per unit 
mass in clusters. § [5] reviews and discusses the papers main findings. 



A flat cosmology with H = 70 km s 1 Mpc 1 , Q 



Lambda 



0.7, and Q 



M 



0.3 is used. 



2. OBSERVATIONS 



Three catalogs are described here, all of which were generated from the four Deep 
fields of the Canada-France-Hawaii Telescope Legacy Survey: the Supernova Legacy Survey 
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catalog of type la supernova^], the Ilbert et al. (2006) photometric redshift galaxy catalog, 
and the Olsen et al. (2007) optical galaxy cluster catalog. 



2.1. The CFHT Supernova Legacy Survey Catalog 

SNLS images four 1 deg 2 Deep fields (D1-D4) every three to four nights (when visible) 
in four MegaCam filters [gu^M^M-, z m) to a depth i M ~ 25. At the end of its five year 
program, the SNLS wil l have discovered ~ 500 type la supernova with its imaging and 



spectroscopic programs (IHowell et al. 20051 ). and be a valuable contribut or to collaborative 



efforts in constrain ing the dark energy equation of state parameter w (lAstier et al. 2006 



Spergel et al. 20071 ). For this project we use the 343 SNela spectroscopically identified prior 



to 2007 September 29. 



2.2. Ilbert Photometric Redshift Galaxy Catalog 

The Ilbert et al. (2006) galaxy photometric redshift catalog^ covers the four SNLS Deep 
fields. They incorporate VIMOS VLT Deep Survey spectroscopic redshifts to calibrate the 
spectral energy distribution (SED) fitting routine, resulting in photometric redshifts with an 
accuracy of (Taz/(i+z) = 0.029 for i AB < 24, and a fractio n of catastrophic errors of 1% at 
17.5 < iab < 21.5, increasing to 10% at 23. 5 < i AB < 24 Jllbert et al. 2006h . In optimizing 



the photo-z calculation, the accuracies of the SEDs types are compromised (Ilbert, private 
communication), and the distribution of SED types is discontinuous. To solve this problem, 
we use 51 SEDs interpolated from Coleman et al. (1980) and Kinney et al. (1996) templates^, 
and fit them to the apparent magnitudes and photometric redshifts of the catalog galaxies. 
We then calculate K corrections and absolute magnitudes, and estimate galaxy stellar masses 
and star formation rates using fits of this library of SEDs to the models of Buzzoni et al. 
(2005), which we find to agree well (within a factor of 2) with those derived from PEGASE 



models (ISullivan et al. 2006bl ). 



To ensure catalog purity, we restrict the catalog as recommended by Ilbert et al. (2006). 
For example, galaxies must be detected in i M , the flag values are used to mask galaxies 
near foreground stars, the number of bands used for redshift fit must be at least 3, and 



1 http: / /legacy.astro.utoronto.ca 

2 http://terapix. iap.fr 

3 http: / / www.astro.uvic.ca/~gwyn/cfhtls 
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galaxies with a second peak in their redshift Probability Distr ibution Function indicate a 
likely catastrophic failure, thus are excluded (jllbert et al. 20061 ). In addition we impose the 
limits iyi < 25 and z < 1.2, and that the "object" parameter must be equal to zero indicating 
the object is a galaxy. 



2.3. Olsen Optical Galaxy Cluster Catalog 

The Olsen et al. (2007) cluster catalog derived from the Terapix data release in August 
2005 (T0002) includes photometric redshifts and an optical grade (A-D). Olsen clusters are 
restricted to those of grade A, meaning a concentration of similarly colored galaxies is visible, 
resulting in 18, 17, 16, and 10 clusters in Dl-4. This optical cluster catalog was chosen over 
those of other selection techniques as it provides consistent completeness for all four SNLS 
Deep field areas. 

Olsen et al. (2007) find their cluster redshifts are slightly overesti mated for z < 0.6 , 



and underestimated for z > 0.7, with a standard deviation of ~ 0.1 (jOlsen et al. 20071 ). 
For the cluster redshift we instead use the peak of the redshift distribution from Ilbert 
catalog E/S0 and Sbc galaxies along the cluster's line of sight (within 40 arcseconds). For 
distributions which plateau instead of peak, the plateau's central redshift is used. To evaluate 
corrected cluster redshift precision we apply the same procedure to the VVDS spectroscopic 
redshift catalog^. For the four Dl clusters with VVDS galaxies along their line of sight, the 
differences between photometric and spectroscopic redshift peaks have a standard deviation 
of &&z/{i+z VVDS ) = 0.023, indicating good agreement. 



2.4. Additional Catalog Processing 

SNe la host galaxies are identified as the closest neighbor in the Ilbert catalog - except 
when the offset difference between the nearest two is less than 2 arcseconds, in which case 
redshift is used to discriminate. 85 SNe la have no neighbor within 5 arcseconds (maximum 
host offset), and cannot be used; this mainly includes SNe la in hosts of %m > 25, and on 
the masked regions covering ~20% of field area. For the SNe la with potential Ilbert catalog 
hosts, iterative outlier rejection is applied to the residual dispersion between host photometric 



http : / / www.oamp . fr / virmos / 
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and SNIa spectroscopic redshiftsj The photometric redshift accuracy when performed for 
all fields combined is <7az/(i+z S n) = 0.030 (consistent with Ilbert et al. 2007). However, if 
outlier rejection is applied to each Deep field separately, we find that each has a different 
photometric redshift uncertainty: cr D i = 0.028, am = 0.024, odz = 0.030, a da = 0.050. In 
this way, 26 SNela are rejected as outliers, dropping the total number to 232 SNela. 

The SN la rate per year, SNR Ia , is calculated for every galaxy based on the two compo- 
nent model, parametrized by SNRj a = A M + B M where M is the stellar mass of the galaxy 
(in solar masses) and M is the star formation rate (in solar masses per year). The most re- 
cent A and B values of Sullivan et al. (2006a) are used: A = 5.3 ± 1.1 x 10~ 14 SNe y" 1 M" 1 
and B = 3.9 ± 0.7 x 10~ 4 SNe y _1 (M y -1 ) -1 . The expected SNIa rate per year for each 
galaxy is corrected to reproduce the total number observed in the SN la sample (excluding 
SNela with no identified host). This correction is performed for separate redshift ranges 
and Deep fields to account for incompleteness in observations; values are given in Table [TJ 
The result is the number of SN la expected in each galaxy over the observing period. These 
calculations are also repeated using the A (mass) component only. 



3. PARAMETRIZED CLUSTER-LIKE ENVIRONMENTS 

To begin we use a continuous clustering strength estimator to compare the SN la rate per 
unit mass in and out of clustered environments. This parameter quantifies the significance 
(£) of finding N F neighbor galaxies the environment of a galaxy or SNIa, given number 
expected from the background distribution. S is calculated for a cylindrical volume envi- 
ronment of diameter D and redshift depth ±er Az/(i+z)(l + z ) (°" from § 12.41) centered on the 
object, as in equation [T] where N F is the number of galaxies in the environment, N F is the 
number of galaxies within ±<7az/(i+z)(1 + z) in the Deep field, Ae = tiD 2 /4: is the aperture 
area, a nd Ap is the area of the field after allowing for the area occulted by foreground star 



masks (lllbert et al. 20061 ). Since N F is computed at the same redshift as the object being 



studied, incompleteness at high redshift is automatically compensated for. 

e _ Ne-N f (A e /A f ) 
^N f (A e /A f ) 

The advantage of this parametrization is that a strict cluster definition is avoided - any 



5 A sigma of three and fractional convergence threshold of 0.01 was used; average number of iterations 
was 4 
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desired scale of clustering can be explored by altering the aperture diameter D. Environ- 
mental significances are calculated using a range of clustering scales (diameters D ranging 
from 0.1 to 1.5 Mpc) both for field galaxies (0.1 < z < 1.1), and also for the 232 SNela 
surviving the cuts in section 12.41 Normalized cumulative significance distributions for both 
field galaxies and SNela are shown in Figure[TJfor D = 0.6 Mpc. A Kolmogorov-Smirnov test 
shows that the two samples are statistically indistinguishable: SN hosts appear to be drawn 
from the same population as field galaxies with respect to clustering in their environment. 

Let us now turn to the most clustered environments, by defining a significance limit 
(£x%) to isolate the top 10%, 5%, and 1% most significant galaxy environments - the high 
significance group (HSG). Figure [1] shows the HSG cutoffs for S 10 % as an example; actually 
the HSG is determined for each diameter D and each Deep field separately, then the HSGs 
for all are fields combined for a given D. We use summed Poisson probabilities to compare 
the number of SNela in the HSG (N b s ) to the number expected (N exp ), which is the sum of 
the expected number of SNe la in each HSG galaxy (from § 12 .4p. The Poisson probability o f 



observing x = N Q b s given fi = N exp is expressed by equation 121 (IBevington fc Robinson 20031 ). 
and the summed Poisson probabilities for x > fi and x < fi are given in equations [3] and H] 
respectively. 

P p (x ;/ i)=/i»i^ (2) 



Psum(x, oo;//) = P p (x';ii) (3) 

x'=x,x+l,... 



=3-1 



p SU m(o,x;v)= p p( x ''^) ( 4 ) 



=0,1, 



As expressed by the morphology- dens ity relation, galaxy clusters are dominated by 



early-type galaxies (jPostman Sz Geller 19841 ). so Poisson probabilities are also calculated for 
the subset of early-type field galaxies and SNela hosts in the HSG. Results for a represen- 
tative sample of environment diameters in Table [2] show the number of HSG SNIa to be 
consistent with expectations of the two-component model over a range of D and Thus 
we conclude this clustering parametrization method does not show an influence of clustering 
on SNIa events. Neither including the SNela outliers rejected in § 12.41 nor increasing the 
environment redshift depth to ±2cr(l + z), affects this conclusion. 
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4. SN la IN OLSEN CATALOG CLUSTERS 

Here we identify SNela and galaxies in grade A clusters from the Olsen et al. (2007) 
cluster catalog. We use Poisson probabilities (§ 14.11) and a direct calculation of SNela rates 
in clusters (§ 14. 2[) to look for an influence of clustering on the SNIa rate per unit mass. As 
the SNLS detection efficiencies of Neill et al. (2006) are valid for 0.2 < z < 0.6, we restrict 
cluster, galaxy, and SN la redshifts to this range to use of them. This decreases the catalogs 
to 30 clusters, 70587 galaxies, and 109 SNela. We note Olsen et al. (2007) and Ilbert et al. 
(2006) use different Terapix data releases (T0002 and T0003 respectively), and the release 
used by Olsen et al. has more masked regions. However the differences in field effective areas 
are < 10%, so likely only < 3 clusters of Ilbert catalog galaxies are missing from the Olsen 
catalog. 

SNe la and galaxy members of clusters are identified as are neighbors in the environment 
of a galaxy in § El except the volume is centered on the cluster coordinates. Since the cluster 
filter used by Olsen et al. (2007) has a profile with core radius r c = O.lSS/i^ 1 Mpc and 
cut off radius of r co = 1.33/tf 5 Mpc, results for 0.4, 0.8, and 1.5 Mpc will be presented as 
representative. For galaxies a redshift depth of ±.2^J{a^ z /(i +ZSN )) 2 + (crAz/(i+z VVDS )) 2 (l+^c) 
is used, a convolution of uncertainties in cluster redshifts from § !2.3l and galaxy redshifts from 
§ 12.41 As SNe la have spectroscopic redshifts, a redshift depth of simply ±2a/± z /(i +ZvvDS ) is 
appropriate. Images in Figure [2] and data in Table [3] present 6 SNe la identified in Olsen 
clusters. While the 3 SNela within D = 0.8 Mpc are probably physically associated with 
the clusters, this cannot be said for the remaining three. The number of field SNela with 
0.2 < zsn < 0.6 predict ~ 1.8 SNela would randomly appear in the regions between r = 0.4 
and r = 0.75 Mpc and ±2o"A z /(i+ ZyvDS ). The probabilities of all 3 being random and all 
3 not being random associates are both < 20%. As it remains likely that at least one is 
physically associated with a cluster we include them in our results, but remind the reader 
that interloping SNela (and galaxies) will add to the uncertainties for D = 1.5 Mpc. 

4.1. Summed Poisson Probabilities 

Summed Poisson probabilities (defined in § [3]) of observing these cluster SNe la are com- 
puted from the number expected in the identified member galaxies, for all galaxy types and 
early-types only. Table H] shows these observations are consistent with the two-component 
SNIa rate model. In fact, > 6 SNela would have to have been observed within D = 0.4 
Mpc of galaxy clusters (> 8 or for D = 0.8 Mpc) for Psum < 0.05. This constrains the 
SNIa rate in clusters to agree with the two-component model to within a factor of two. 
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Obviously, to detect an effect on SNela rates due to clustering (over and above the 
expectations of the two component model and morphology-density relation) lurking in this 
data, we should at least detect the two component model. So can we rule out the single 
component model (mass or "A" component only) from these cluster observations? The 
answer appears to be no. Table H] shows that the number of observed cluster SNIa is actually 
consistent with bot h models. Surveys to test for the A+B model or more complicated delay 
time distributions (IMannucci et al. 20061 ) will require a larger survey (two to three times the 
area or duration); unfortunately, the final SNLS data set will not be adequate for this. 



4.2. SNIa Rates in Clusters 

Here we use detection efficiencies from Monte Carlo simulations of the SNLS SNIa 
identification pipeline from Neill et al. (2006) to directly calculate the SNIa rate per unit 
mass in clusters, and compare it to that for the field. The corrected number of SN la which 
exploded in a given field in a year, N corr j a , is extracted from equation 3 of Neill et al. 
(2006). As shown in equation [5], Nj a is the total number of supernova observed in a given 
field; S is the number of observing seasons; e yr is the detection efficiency per year (the 
probability of a supernova being detected, sent for spectroscopy, and identified as a type la); 
and Cspec accounts for the fraction of SNela for which spectra are obtained yet remain 
unidentified. [1 + {z)y\ corrects for time dilation at the volume- weighted average redshift of 
the survey, and (z)y = 0.46. Detection efficiencies from Neill et al. (2006) and the number 
of observing seasons for each field up to January 2007 are in Table [51 Approximations to 
Poisson uncertainty from Gehrels (1986) determine the upper and lower limits on Nj a at 
the 0.84 confidence level (corresponding to la), which are substituted into equation [5] to 
determine AN corr j a . 

N corr , Ia = N ' a/S [l + (z) v ] (5) 

tyr^SPEC 

Table [6] contains the resulting corrected SN la rate per unit mass per year in clusters 
and the field. With only 3 identified cluster SNela, statistical uncertainties dominate this 
calculation; systematics, mainly the error in galaxy mass calculations, are likely another 
~30%. Thus, our SNIa rat e is consistent with b oth the rate in early-type galaxies, 5.3 ± 



1.1 x 10 SNe y M ( ISullivan et al. 2006aT). and the low redshift cluster rate from 



WOOTS, 9.8±|;| ± 0.9 x 10" 14 SNe y" 1 M" 1 flSharon et al. 2007h . 



There are two factors not considered which could affect SNIa detection efficiencies in 



cluster galaxies. First, SN la detection efficiencies decrease in brighter hosts (INeill et al. 20061 ) 
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and the brightest galaxies are early-type. Second, the SNIa detection efficiency decreases 
for fainter, lower stret ch SN la, and these faint SN la occur preferentially in early-type hosts 
( ISullivan et al. 2006al ). Since cluster galaxies are predominantly early-type, both of these 
effects dominate in clusters: the first effect decreases the number expected in clusters by 
~ 15%, but quantifying the second would require more detailed completeness simulations. 
Both effects would cause us to underestimate the rate of SNela in clusters relative to the 
field. 



4.3. SNIa Rates in E/SO Cluster Galaxies 

To avoid these effects and the morphology-density relation, we limit the galaxy sample 
to two subsets: all early-type galaxies, and the brightest population of early-type galaxies 
(those with My < —23.0 like brightest cluster galaxies, BCGs). This has the added benefit 
of rejecting interlopers misidentified as cluster members. Two cluster SNe la have early-type 
hosts, and the host of 03Dlax is brighter than My = —23.0. Detection efficiency corrections 
are performed as described above, with the final results and rates presented in Table [6j 

Although these samples are more sensitive to the two detection efficiency biases affecting 
early-type galaxies (§ 14.21) . limiting all galaxies to early-types minimizes differences between 
the clusters and the field and results in a more meaningful test. The results, although not 
statistically significant, are suggestive of the rate enhancement in cluster over field early-type 
galaxies established by Mannucci et al. (2007). An enhancement in BCG-like galaxies would 
be consistent with the findings of an enhanced SNIa rate in radio loud elliptical galaxies 



(IDella Valle et al. 20051 ). as these are usually the brightest cluster galaxies. Future work 
will investigate the rates of SNLS SNela in radio galaxies using existing radio catalogs for 
the Deep fields. 



4.4. Effects of Altering Data Constraints 

Including grade B clusters from the Olsen et al. (2006) catalog increases the total 
number of clusters to 65, with one new cluster SNIa identified (03Dlfb at z$n = 0.498 in 
an E/S0 host of My = —21.9 with cluster offset 57.8"). This does not affect the conclusions 
of the Poisson probability experiment (§ 14.11) . The resulting cluster rates for D = 0.8 Mpc 
for all galaxy types, early-types, and the brightest early-types only are 6.7l 3 ; , 6.7^3 1 4 °, and 
6.2lg 9 { 2 xlO -14 SNe y -1 Mq 1 ; slightly below yet well within uncertainties of results with A 
grade clusters only. 
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None of the SNela rejected as outliers in § 12 .41 are associated with Olsen catalog clusters. 
Including them does not affect the conclusions of the Poisson probability experiment, and 
only increases the SNela field rate for all galaxy types to 15.9^ J xlO~ 14 SNe y _1 Mq 1 . 
Including in the experiment all SNela for which no host was identified does not yield any 
new cluster SNe la either. 

Reducing the SNela sample to those included in Neill et al. (2007) avoids a possible 
over-correction for detection efficiencies which might arise from extending the sample to 
later times, if the survey became more efficient (we estimate any improvement in detection 
efficiency to be small). This decreases the number of SNela to 40, with one cluster SNela 
(03Dlax). This does not alter the overall conclusions of these experiments. 

Lastly, the inclusion of SNLS SNela candidates which did not receive spectroscopic 
confirmation (past maximum light, or detected as field season was ending), but which were 
given photometric types and redshifts using the techniques of Sullivan et al. (2006b), does 
not yield any new cluster SNe la or affect the conclusions of these experiments. 

5. CONCLUSIONS 

Type la supernova, galaxy, and cluster catalogs generated from the first four years of 
CFHTLS Deep survey data were combined to search for an influence of clustering on the 
SNIa rate per unit mass. To avoid cluster-specific detection efficiencies and the inclusion of 
interloping galaxies as cluster members, we also considered subsets of regular and BCG-like 
early-type galaxies only. Results are dominated by the statistical uncertainties in identifying 
only three probable and three possible cluster SNela, and consistent with the results of 
low redshift cluster SNela rate studies. To capitalize on SNLS's inclusion of both field and 
cluster environments, we used the continuous clustering strength parameter "significance" 
to compare the SN la rate per unit mass in and out of clustered environments, but did not 
find evidence of a clustering influence on the SN la rate per unit mass. Future work includes 
incorporating existing radio and infrared catalogs to investigate the rates of SNLS SNela 
in radio loud and infrared luminous galaxies, and the VIRMOS spectroscopic redshifts for 
Deep field galaxies to explore the SNe la rate in small groups. 
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and Henry McCracken for early access to and correspondence regarding their photometric 
redshift galaxy catalog, Lisbeth Olsen for early access the optical cluster catalog, Jon Willis 
and Dan Maoz for helpful conversations, and our anonymous referee for constructive cor- 
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Table 1: SNIa rate correction factors. 



Redshift Range 


Dl 


D2 


D3 


D4 


0.0-0.5 


1.431 


1.016 


1.538 


0.896 


0.5-0.7 


1.203 


0.511 


0.940 


0.999 


0.7-0.9 


0.523 


0.194 


0.418 


0.616 


0.9-1.1 


0.078 


0.065 


0.140 


0.118 


1.1-1.3 


0.000 


0.000 


0.000 


0.000 
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Table 2: Results of the clustering parametrization method for Dl-4. 



Diameter 




Sio% 






£5% 






£1% 




(Mpc) 


N obs 


N 

1 > exp 


PsUM 


N obs 


N 

1 ' exp 


PsUM 


Nobs 


N 

1 ' exp 


PsUM 


0.4 


21 


26.9 


0.15 


11 


15.0 


0.18 


3 


3.98 


0.44 


0.6 


27 


26.5 


0.48 


16 


14.7 


0.41 


5 


3.41 


0.26 


0.8 


21 


25.3 


0.23 


12 


13.8 


0.38 


5 


3.15 


0.21 


1.0 


23 


24.8 


0.41 


10 


13.2 


0.23 


6 


3.13 


0.10 


1.5 


21 


22.3 


0.44 


11 


12.0 


0.46 


2 


2.98 


0.43 


For early-type galaxies 


and SN la hosts 


only: 










0.4 


5 


7.43 


0.25 


3 


3.77 


0.48 





0.72 


0.49 


0.6 


6 


6.86 


0.47 


3 


3.27 


0.59 





0.60 


0.55 


0.8 


5 


6.33 


0.39 


3 


3.08 


0.63 


1 


0.49 


0.39 


1.0 


3 


5.98 


0.15 


3 


2.86 


0.54 


1 


0.59 


0.45 


1.5 


2 


5.72 


0.08 


2 


2.76 


0.48 





0.61 


0.54 
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Table 3: Cluster SNela details. 



SNIa 


SNIa 


SNIa 


SNIa 


SNIa 


host 


host 


host 


cluster 


cluster 


SNLS ID 


RA 


Dec. 


%spec 


stretch 


type 


My 


offset 


Zphot 


offset 


03Dlax 


02 h 24 m 23 s .32 


-04°43'14".41 


0.496 




ES/0 


-23.64 


1.97" 


0.53 


10.8" 


06Dlkg 


02 h 24 m 32 s .57 


-04°15'02".0 


0.323 


1.21 


Sbc 


-20.31 


1.69" 


0.30 


47.2" 


05Dlby 


02 h 24 m 35 s .45 


-04°12'04".2 


0.299 


0.99 


Sbc 


-21.17 


0.86" 


0.30 


138.6" 


04Dlpg 


02 h 27 m 04M6 


-04°10'3i".4 


0.515 


1.05 


Sbc 


-19.62 


0.18" 


0.51 


85.2" 


05D3mq 


14 h 19 m 00 s .40 


+52°23'06".81 


0.246 


0.90 


ES/0 


-21.85 


4.94" 


0.25 


9.99" 


07D3af 


14 h 19 m 05 s .01 


+53°06'08".98 


0.356 


0.98 


Scd 


-18.61 


0.35" 


0.30 


128.7" 
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Table 4: Summed Poisson probabilities for cluster SNela for Dl-4. 



Cluster A+B Components 


A Component 


Diameter N obs N exp 


PsUM 


N 

1 ' exp 


PsUM 


0.4 2 2.38 


0.57 


3.21 


0.38 


0.8 3 3.51 


0.54 


4.53 


0.34 


1.5 6 5.97 


0.55 


7.30 


0.41 


For early-type galaxies only: 






0.4 2 1.97 


0.59 


2.79 


0.47 


0.8 2 2.60 


0.52 


3.67 


0.29 


1.5 2 3.69 


0.29 


5.23 


0.11 
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Table 5: Parameter values for equation [5j 



Field 


Seasons 


a 


CsPEC 


Dl 


4.0 


0.3 


0.94 


D2 


4.0 


0.22 


0.88 


D3 


4.0 


0.31 


0.80 


D4 


4.1 


0.31 


0.69 



From Neill et al. (2006). 
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6: Corrected SNela 


rates 


per unit mass for Dl-4. 




Galaxy 


N Ia 


N T 

1 * corr,la 


Stellar Mass 


SNIa Rate 


Set 




(SNe y" 1 ) 


(xlO 14 M Q ) 


(xl0~ 14 SNe y" 1 M" 1 ) 


Clusters, 0.4 Mpc 


2 




0.33 


o q+19.4 

O.0_ 4 g 


Clusters, 0.8 Mpc 


3 




0.48 


8-5™ 


Clusters, 1.5 Mpc 


6 




0.78 


10.4lft 


Whole Field 


109 


i69.il":! 


11.4 


14.811:1 


For early-type galaxies only: 






Clusters, 0.4 Mpc 


2 


2.81ft 


0.29 


9-6™ 


Clusters, 0.8 Mpc 


2 


2.81ft 


0.38 


7 +16.8 
I -^-4.2 


Clusters, 1.5 Mpc 


2 


2.81ft 


0.56 


5.0™ 


Whole Field 


23 


35.1™ 


6.23 


5.6111 


For the brightest early-type galaxies 


only: 




Clusters, 0.4 Mpc 


1 




0.13 


9.8™ 


Clusters, 0.8 Mpc 


1 


1 3+ 6 ' 1 


0.15 


8.7™ 


Clusters, 1.5 Mpc 


1 


1 3+ 6 - 1 


0.19 


7.01| 3 8 ° 


Whole Field 


2 


2.6±f : ? 


1.02 


2.5^:| 
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X 
CO 
CO 

CD 



o 



u 

CO 




All Fields Combined 
Diameter 0.6 Mpc 

Galaxies 

SNIa 

KS D 0.0614611 

KS P 0.335906 
_. i I i ._ 



2 4 
Significance 



Fig. 1. — Cumulative significance distribution for galaxies and SNela, with environment 
diameter 0.6 Mpc, for Dl-4. Kolmogorov-Smirnov maximum difference (KS D) and prob- 
ability (KS P) show the null hypothesis cannot be ruled out. Dashed lines represent £io%- 
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Fig. 2. — SNLS SNela identified in galaxy clusters from the Olsen et al. (2007) catalog. 
SNela coordinates marked by crosses, clusters represented by circles of diameter 0.4 Mpc (A, 
B), 0.8 Mpc (inner circle of C), and 1.5 Mpc (outer circle of C, D, and E). Images created 
from SNLS 2004 reference images and Skycat (image quality degraded for submission to 
astro-ph) . 



